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Abstract 

Information on su(N) tensor product multiplicities is neatly encoded in Berenstein-Zelevinsky 
triangles. Here we study a generalisation of these triangles by allowing negative as well as 
non-negative integer entries. For a fixed triple product of weights, these generalised Berenstein- 
Zelevinsky triangles span a lattice in which one may move by adding integer linear combinations 
of so-called virtual triangles. Inequalities satisfied by the coefficients of the virtual triangles 
describe a polytope. The tensor product multiplicities may be computed as the number of 
integer points in this convex polytope. As our main result, we present an explicit formula for 
this discretised volume as a multiple sum. As an application, we also address the problem of 
determining when a tensor product multiplicity is non- vanishing. The solution is represented 
by a set of inequalities in the Dynkin labels. We also allude to the question of when a tensor 
product multiplicity is greater than a given non- negative integer. 
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1 Introduction 



The decomposition of tensor products of modules of simple Lie algebras has been studied for 
a long time now. Many elegant results have been found for the multiplicities of the decompo- 
sitions, the so-called tensor product multiplicities. The relatively recent Berenstein-Zelevinsky 
method of triangles [PJ is an example. Although it is a powerful, symmetric method, it is not 
explicit: triangles are constructed according to certain rules, and their number is the required 
tensor product multiplicity. Here we show that a generalisation of these Berenstein-Zelevinsky 
(BZ) triangles allows us to work out a very explicit expression for the multiplicities. A tensor 
product multiplicity is expressed as a multiple sum, counting the number of integer points in 
a particular convex polytope, to be defined below. BZ triangles and our results pertain to the 
A-series; A r = su(r + 1). We will sometimes write su(N) with N = r + 1. 

We are interested in describing decompositions of tensor products of irreducible highest 
weight modules of simple Lie algebras. They are usually written 

M A <g> M„ = T x /M v , (1) 

where M\ is the module of highest weight A. T\ tfl u is the tensor product multiplicity. We shall 
study the equivalent but more symmetric problem of determining the multiplicity of the singlet 
in the expansion of the triple product 

M\ <S> (g) M v D T X ^M . (2) 

If v + denotes the highest weight conjugate to v, we have T\ ifJltV = T\^ + . We will use the 
shorthand notation A <8> (jl ® v to represent the left hand side of (Q) . 

An su(3) BZ triangle, describing a particular coupling (to the singlet) associated to the 
triple product A ® fx ® u, is a triangular arrangement of 9 non-negative integers: 



m 13 
"12 l 2 3 
m 2 3 mi 2 

"13 l\2 "-23 ^13 



(3) 



These integers are related to the Dynkin labels of the three integrable highest weights by 

mi3 + n 12 = Ai , ni3 + lyx = Ml , ^13 + ^12 = ^1 , ,^ 

"T-23 + "-13 = A 2 , n 2 3 + ^13 = ^2 , ^23 + ™43 = ^2 • 

We call these relations outer constraints. The entries further satisfy the so-called hexagon 
conditions 

"12 + f"23 = "23 + mi2 , 

mi2 + ^23 = "T-23 + ll2 , (5) 
^12 + "23 = ^23 + "12 

of which only two are independent. They say that the length of opposite sides of the hexagon 
must be equal, if the length of a segment is defined to be the sum of the two integers associated 
to its endpoints. An su(3) BZ triangle is thus composed of one hexagon and three corner points. 
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For su(4) the BZ triangle is denned in a similar way, in terms of 18 non-negative integers: 



mu 
"12 hi 
m 24 m 13 

"•13 ^23 "-23 hi 

m 34 m 23 mi2 

ni4 l\2 " 2 4 ^13 "-34 ^14 



(6) 



related to the Dynkin labels by 

m\A + "12 = Ai , nu + Z12 = Ml > 'i4 + "H2 = v\ , 

"1-24 + "-13 = A 2 , n 2 4 + hz = M2 , ^24 + ™13 = ^2 , (7) 

m 34 + nu = A3 , "34 + Z14 = yu 3 , Z34 + m 14 = z^ 3 . 
Furthermore, the sit(4) BZ triangle contains three hexagons: 

"12 + m 2 A = m 13 + n 23 , n i3 + Z 23 = Z 12 + n 24 , hi + "23 = l\3 + "34 , 

"12 + ^34 = ^23 + "23 , "13 + "T-34 = "24 + "*23 , "23 + "1-23 = "H2 + "34 , (8) 

"i-24 + h3 = hi + mi 3 , m 34 + l 12 = I23 + m 23 , I13 + ™23 = hi + m 12 . 

It is a general feature for any N that only two out of the three hexagon identities associated to 
a particular hexagon are independent. 

The su(N) generalisation is obvious; the triangle is built out of (JV — 1)(N — 2)/2 hexagons 
and three corner points. Simple examples of lower rank BZ triangles and their applications may 
be found in Ref. . 



2 Generalised and virtual Berenstein-Zelevinsky triangles 

The generalisation of the BZ triangles we shall consider is obtained by weakening the constraint 
that all entries are non-negative integers to arbitrary integers, negative as well as non-negative. 
The hexagon identities and the outer constraints are still enforced. A triangle will be called a 
true BZ triangle if all its entries are non-negative. 

We consider a generalised BZ triangle associated to su(r + 1). Denoting the number of 
entries E r and the number of hexagons H r , we have 

E r = \r{r + 1) , H r = K(r - 1) . (9) 

For a given triple product A ® /j, v, the set of associated triangles spans an .ff r -dimensional 
lattice. Each hexagon corresponds to two independent constraints on the triangle entries while 
there are 3r outer constraints. This leaves 

E r - (2H r + 3r) = H r (10) 

parameters labelling the possible triangles. Among these, only a finite number are true BZ 
triangles. This number is precisely the tensor product multiplicity of the triple coupling. For 
example, when the singlet does not occur in the decomposition of the triple product, there are 
no true BZ triangles in the lattice. 
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A special class of generalised BZ triangles is associated to the triple product 000(8)0. We say 
they have weight (A, fj,, v) = (0, 0, 0). According to the general argument above, H r = \r(r — 1) 
such triangles are linearly independent. We shall call them virtual triangles, and denote them 
using V. It is natural to exclude the triangle with all entries equal to zero from the set of 
virtual triangles. It is the unique true BZ triangle in the lattice associated to the triple product 
(g> <X> 0. In the cases of su(3) and sii(4), virtual triangles have appeared in Ref. || and Ref. 
Q, respectively. 

A convenient basis for the virtual triangles is given by associating a simple distribution of 
plus and minus ones (written 1 and 1 = — 1, respectively) to a given hexagon. All other entries 
are zero. The distribution is 

1 I 11 

1 1 (11) 

1111 
1 

Thus, a basis virtual triangle will always have 6 entries equalling —1, and between 3 and 6 
entries equalling +1. The number of +1 entries depends on where the associated hexagon is 
situated in the generalised BZ triangle and on the rank of su(r + 1) (3 entries equal +1 for 
su(3) only). That these virtual triangles are indeed linearly independent is obvious. 

There are no virtual triangles in the case of su(2). In the case of su(3) there is one basis 
virtual triangle 

1 

V= j 1 1 1 (12) 
1111 

while in the case of su(4) the three basis virtual triangles Vi, V2 and V3 are 

1 
11 

Vl= 1 - - l V 2 = - 1 - ° 

1 1111 2 1110 

10 110 

000000 11I100 





V 3 = „% j 1 ! (13) 

1 1 

1111 

The generalisation to higher rank su{r + 1) is straightforward. In Section 3 we shall use another 
choice of indices on V. 

We are now in a position to generate all generalised BZ triangles associated to a given triple 
coupling. Once a single generalised triangle has been found, the lattice of triangles associated 
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to the triple coupling is spanned by adding integer linear combinations of the virtual triangles. 
The choice of initial triangle is not important. We shall denote the integer coefficients linear 
coefficients. We emphasise that the euclidean spaces spanned by the lattices are all of dimension 
H r , i.e., the dimension is independent of the triple coupling and depends only on the rank of 
su(r + 1). 

Negative entries in (generalised) BZ triangles appear also in Ref. ||. That work is mainly 
devoted to the construction of tensor-product generating functions. A new method is proposed 
based on elementary solutions to certain sets of linear equations related to the BZ triangles, and 
for sit(3) one of these solutions corresponds to a triangle with negative entries. The appearance 
of negative entries is expected to be a general feature for higher su(N) as well. The elementary 
solutions are closely related to the so-called elementary couplings. 



3 Polytopes, multiple sums and tensor product multiplicities 

We shall now focus on the ff r -dimensional linear coefficient space and seek an algebraic descrip- 
tion of the tensor product multiplicities. The latter are computed by counting true BZ triangles. 
Demanding that the entries of a true BZ triangle be non-negative, we obtain E r inequalities 
the linear coefficients must satisfy. The inequalities depend on the choice of initial triangle, 
and they correspond to a polyhedral combinatorial expression for the multiplicities in linear 
coefficient space. The structure of the basis virtual triangles (cf. (pT|), (12) and (|lq)) ensures 



that all linear coefficients have upper as well as lower bounds. The polyhedron is therefore 
bounded and such a polyhedron is called a polytope. It is easily seen to be convex. 

In order to specify the polytope, we must find an initial triangle of weight (A, (J,, v). It 
is convenient to break the symmetry among the 3 weights, and first look at the unique true 
triangle of weight (A, fj,, A + + fi + ). For su(3), this triangle is 



Ai 
/ii 
A2 A2 

m H2 



(14) 



and the generalisation to su(r + 1) is clear. Every highest weight v in a coupling A ® [i ® v 
satisfies 

r 

v = A + + - ^2 nidi , (15) 
i=l 

with Hi € where on is the i-th simple root. The coefficients are conveniently expressed 
using dual Dynkin labels. A weight A can be written 

A = ^>A 4 = j^Xat , (16) 

i=l i=l 

where {A*} and {ctf} are the sets of fundamental weights and simple co-roots, respectively. The 
A* are the dual Dynkin labels, while the ordinary Dynkin labels are the Aj. For simply-laced 
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algebras, like su(N), a% is identical to the co-root a( (with standard normalisation a 2 = 2, for 
a a long root). Taking the scalar product of (15) with A* therefore gives 

n, = (A+y + Gu+y-i/ . (17) 

Generalised triangles of weight (0,0, en) are also easily constructed. An su(3) example is 



1 

1 1 

I 




(18) 



of weight (0,0,02)- So, one can find a generalised triangle of weight (A, fj,, u) by subtracting 
non-negative integer multiples of triangles of weight (0, 0, Oj), such as fllSj) , from a triangle like 
(0), of weight (A, n, A+ + 

The result for su(r + 1) is the following generalised BZ triangle associated to the triple 
product A <g> fj, (g) v: 

N' r 
n r N r 
A 2 K_ t 
hi n r _i N r _i 
A3 A 3 N>_ 2 

/J,l H2 n r -2 A r r-2 







A r -2 













Ml 




n 3 iV 3 






-1 




A r _i 


A r _i ATj 







Ml 





M2 


^ r _ 2 "2 


iV 2 






A/- 


A r 


Ay Ay 


N[ 


Ml 





M2 


/i 3 


Hr-2 /i r _i 


Til Ni 



(19) 

The entries rii, Ni and TV/ are defined by 

rii = X r ~ i+1 + H r ~ i+1 ~ , 

Ni = (1 - Sii)m-i - Hi + Hr-i+l 

= _ A r-i+l + (1 _ - (1 - ^/Z^ + - (1 - fo^" 1 + iS , 

Ni = Vi-Ni 

= \ r - l+l - (1 - 8a)X r - i+2 + (1 - - H r ~ i+1 + f* - (1 - <VK +1 • (20) 

In order to be able to describe the polytope explicitly, we need to label the virtual triangles. 
Our choice is to write them as Vj j or V; depending on where the associated hexagons are 
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situated. The corresponding linear coefficients are denoted djj and rjf. 

* * * * 
■k d r -2,l * ri r -2 * 
★ ★ ★ * ★ ★ 

* d 2) i * ... * r?2 * 

* * * * * * * * 

* dj i * di 2 * * di r _2 * ?7i * 

* * * * * * ... * ★ * * * * 

Here a * indicates an unspecified entry while dj and r/j are the linear coefficients of the virtual 
triangles. They are depicted at the centres of the hexagons associated to the corresponding 
virtual triangles. We have chosen two different notations for the virtual triangles (and their 
associated linear coefficients) to reflect the positions of the corresponding hexagons in the asym- 
metric initial triangle (|l9|). 

Now, denoting the initial triangle To, any triangle in the lattice of general triangles of weight 
(A, fi, v) may be written as 

r— 1 i+j=r— 1 

T = T + E VM + E ''':i V ':i ■ ( 22 ) 

1=1 i,j>l 

The associated polytope of interest is in the i? r -dimensional space spanned by dj and r/i. It 
is bounded by the inequalities requiring that all entries in T are non-negative (whereby T is 
ensured to be a pure BZ triangle). Hence, the position of the polytope depends on the initial 
triangle Tq. Nevertheless, the volume of the polytope, the number of integer points bounded 
by the polytope, is independent of Tq. By construction, this number is the tensor product 
multiplicity T\^ v of the triple coupling A <g) \i ® v to the singlet. 

Using the explicit choice of initial triangle fll9| ) and the basis of virtual triangles discussed 
above, it is simple to write down the inequalities defining the polytope. To illustrate, we list 
the three inequalities given by the three entries located furthest to the right on the bottom line 
of T: 

^ r _i + di, r _ 2 - m > , m - 771 > , Nx + 771 > . (23) 

A similar polyhedral combinatorial expression is discussed in Ref. ||. Convex polytopes con- 
structed there lie in the space of Gelfand-Tsetlin patterns (see e.g. [pL 0]), while ours lie in spaces 
associated to BZ triangles. Hence, for su(N) their polytopes are embedded in the euclidean 
vector space while ours may be embedded in the smaller space 

| R (JV-2)(JV-l)/2_ A 

more universal method of constructing polyhedral combinatorial expressions for tensor product 
multiplicities, generalising that of Ref. || and making sense for any simple Lie algebra, may 
be found in Refs. M. 
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3.1 Explicit multiple sum formula 

As already stated, our polyhedral expression differs from the one discussed in Refs. ]|, 0. 
Its structure allows us to extract an explicit multiple sum formula counting the integer points 
bounded by the polytope. The multiple sum is over the linear coefficients, so different orders 
of summation give a total of H r \ possible representations of the polytope volume. For practical 
purposes, however, there are considerably fewer appropriate summation orders. Let us illustrate 
our procedure for choosing an appropriate order of summation by considering the following 
simple planar example. 

Let a planar polytope be defined by the set of inequalities 

l<x<4, 8 < x + y < 14 , 

6 < y , 4<y-x<8. 1 ' 

The volume or area A of the polytope (the number of integer points bounded by the inequalities) 
can be written in two ways: 

4 min{a:+8,14— x} n min{4,j/— 4,14— y} 

A = 16 = £ £ 1 = £ E 1 • ( 2 5) 

x=l j/=max{6,z+4,8— x} J/=6 z=max{l,j/— 8,8— y} 

The second expression is slightly more difficult to write, since the upper limit 11 on y must 
be calculated from the intersection of faces (lines). Here the bounding lines x + y = 14 and 
y — x = 8 intersect at the point (x, y) = (3, 11). This is a complication to avoid when writing our 
formula, since it will involve many sums. In the first expression the explicitly written lower limit 
y = 6 is redundant since the remaining two intersect at the point (2,6). However, including 
redundant limits does not change the result, so we may choose to keep the limit y = 6. 

It is clear that an important difference between the two orders of summation is that we have 
1 < x < 4, but only 6 < y. For example, if (|24j) is supplemented with y < 15, then the upper 
limit 11 can be replaced simply by 15, and the formula is still valid. While it is true that the 
new formula contains single sums with lower limits greater than upper limits, this is a relatively 
small inconvenience. Those sums simply contribute 0. 

So, when choosing an appropriate order of summation (over the summation variables r\i and 
dij), it is crucial that for any summation variable, all subsequent summation variables have 
upper as well as lower bounds parallel to (or independent of) the one under consideration. 

This is a non-trivial consideration: the procedure does not apply to all polytopes. A trivial 



example is provided by (24) with the inequality x < 4 removed. 

Fortunately, the simplifying procedure applies to our polytope. A simple inspection reveals 
that not all summation orders are appropriate, however. Nevertheless, in the general case and 
in accordance with our procedure, we may express the volume of the polytope as 

*w = (e) (ee) - ( E - E Ue -e) i <m 

where the summation variables are bounded according to 

max{-JV 1 ,di r _ 2 , -A^ + m, ~Mr-2 + di,r-2 - d 2 , r -3 + m} 
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< 77i < minjni,^.! + di jr _ 2 , A r - d ljr _ 2 + V2,n 2 + di, r -2 - V2,N[,N 2 + r/ 2 } , 
max{d/_i ir _; - di-i tT -i-i + dz,r-«-i, -N/'+i + ?7h-i> 

— Mr-i-l + ^Z+l _ (1 _ ^,r-2)^+l,r-Z-2 + ^/,r-«-l} 

< ^ < min{A r _ m - c^ r _z_i + dj-i.r-j + r/ m , 

+ ^,r-«-i - ^+1,^+1 + ?7«+i} , for 2 < / < r - 2 , 
max{cf r „2,i, -^r} < r/r-l < min{A 2 + (f r _ 2j i, ra r , iV r } , 
max{dij_i, -fij-i + dij-i + d 2 j-i - (1 - ^,2)^-2} 

< dlj < min{/ij + dij-i, A r — dij-i + d 2 ,j-i} , for 2 < j < r — 2 , 
max{djj_i + — dj-ij-i, —fij-i + dt,j-i + dj+ij-i — (1 — ^7,2)^+1^-2} 

< d itj < Xr-i+i - djj-i + di+ij-i + , for 2 < i, j, i + j < r - 1 , 
< < A r _ i+ i + , for 2 < i < r - 2 , 

< < min{/ii, A r } . (27) 



From (15) and ( |17| ) it follows that the weights are subject to the condition 

A* + + G ^> , i = l,...,r , (28) 
ensuring the integer nature of the entries, and thus also of the summation limits (p7|). The 



multiple sum formula (26) is our main new result. We now demonstrate its usefulness by 
considering an application and working out a few examples. 



4 An application 

It is of interest to know whether or not a coupling of a certain weight (A, v) exists, without 
having to work out the tensor product multiplicity. Based on our multiple sum formula (p6|) 
and (^7j) one may derive a set of inequalities in the dual and ordinary Dynkin labels of the 
three weights, determining when the associated tensor product multiplicity is non- vanishing. 
To illustrate the method, we discuss the inequalities for su(3) and su(4) and outline their 
derivation. In principle, it is possible to repeat the procedure for higher rank, but even for 
su(4) the derivation is very cumbersome. We believe that similar results exist for all simple Lie 
algebras, and hope to report more general results later. 

The dimension of the linear coefficient space for su(3) is one, so the tensor product multi- 
plicity may be represented by a single sum: 

min{/ii, A 2 ,A 2 +/i 2 -v 1 ,A 1 +^ 1 -v 2 ,-\ 1 +\ 2 +fi 1 -v 1 +v 2 ,X 2 +fi 1 -fi 2 +u 1 -u 2 } 

E !• ( 29 ) 

r 1 =ma,x{0,\ 2 +fj, 1 ~n 2 -v 1 ,-\ 1 +\ 2 +fi 1 ~v 2 } 

The weights are subject to the integer constraint (^). Note that the summation limits are not 
symmetric in the weights. This is simply because we have chosen an asymmetric initial triangle. 
The summation ((2^) is non-vanishing if and only if the upper limit is greater than or equal to 
the lower limit. This condition yields 6 • 3 = 18 inequalities: 

< \i,fj,i,Vi , for i = 1,2 , 
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maxIA 1 - A 2 + // - -A 1 + A 2 - //} < v x < A 2 + /i 2 , 

max{-A 1 + A 2 - / u 1 + ^ 2 ,A 1 - / u 2 ,-A 2 + ^ 1 } < v 2 < A 1 + / u 1 , 

max{-A 2 - /i 1 + ,u 2 , -A 1 + A 2 - ^ 2 } < v 1 - v 2 

< minfA 1 -^ 1 +/i 2 ,-A 1 + A 2 + /i 1 } .(30) 

Expressing the inequalities in terms of ordinary Dynkin labels, one should bear in mind 
that the summation variable increases in steps of one while the quadratic-form matrix involves 
a factor of 1/3. A similar factor 1/N is present for higher rank su(N). 

In the case of sit(4) the tensor product multiplicity may be written as a triple sum: 

min{/ii,A 3 } min{A 2 +d,A 1 +/t 1 -i/ 3 ,-A 1 +A 2 +^ 1 -!/ 2 +^ 3 } 

d=0 »?2=max{(i,-A 1 +A 2 +/t 1 -^ 3 } 
min{A 3 +^ 3 -i/ 1 ,^ 2 +rf,A 3 -c(+»y2,A 2 +/i 2 -^ 2 +c(-»y2,A 3 +^ 2 -^ 3 +^ 1 -!/ 2 ,-A 2 +A 3 -^ 1 +^ 2 -!/ 1 +f 2 +»)2} 

E i. 

7/i=max{A 3 +^ 2 — /i 3 — f 1 ,^,— A 2 +A 3 — /i 1 +/i 2 — z^ 2 +j^ 3 +t?2 5 — ^i+<^+»?2} 

(31) 

The weights are subject to the integer constraint (|2g|). For a multiple sum like this, the in- 
equalities are obtained by first considering the interior summation over 771 , leading to 6 • 4 = 24 
inequalities which may depend on the remaining two summation variables: 4 of them depend 
only on the weights, 6 depend on d but not 772, while 14 depend on r? 2 . Treating the latter 
in the same way as the upper and lower bounds on the r/ 2 -summation, we obtain a total of 
13 • 6 = 78 inequalities from the ^-consideration. Repeating the procedure for d leads to a 
total of 54-10 = 540 inequalities in addition to the ones already derived from the rjx- and 
^-considerations. This huge set of inequalities may be reduced considerably and we find the 
following constraints on the Dynkin labels (expressed in terms of dual as well as ordinary ones) 



< Aj, fj,i, Vi , for i = 


1,2,3 


5 






max{A 3 - n 1 , A 3 - A 3 - fi 1 + m, -A 1 + /j 3 , -A 1 + A x + /x 3 


-^3} 


< 


v 1 


< A 3 + ii 3 


max{|A 2 - n 2 \, |A 2 - A 2 - [J 2 + (Jq\, |A x - A 3 + /1 1 - 


V!} 


< 


2 

V 


< \ 2 + fj 2 


max{A x - //, A 1 - Ai - fJ? + p 3 , -A 3 + /j 1 , -A 3 + A 3 + /j 1 


-Ml} 


< 


V 3 


< A 1 + m 1 


max{A 2 - A 3 - //, -A 1 + fi 2 - /j 3 , A 1 - A 2 + fj} 




< 


v x 


- fi 


< min{A 2 - A 3 + //, A 3 + fi 2 










max-f-A 1 + A 3 + [i 2 - /J 2 , A 2 - A 2 - // + ijl 3 , A 1 - A 3 + /i 2 


-M 2 , 








A 2 - A 2 + /j 1 - fi 3 , -A 2 + ix 2 - 112, A 2 - A 2 


-A 


< 


v 2 


- f 2 


< min{A 2 - A 2 + /j 2 , A 2 + /j 2 


-M 








max{A 2 - A 1 - /x 3 , -A 3 + i? - /j, 1 , A 3 - A 2 + fi 3 


-m 2 } 


< 


v 3 


- ^3 


<min{A 2 -A 1 + /i 1 ,A 1 + / Li 2 


-m 1 } 








max{-A x + A 3 - /i 2 , -A 2 - // + /U 3 , A 2 - A 2 + ^ 2 




< 


v 1 




< min{-A 1 + A 3 + /j 2 , A 2 - // + /u 3 , A 2 - A 2 + 1? 


- M 









A recent discussion Q includes a brief review of the problem of determining when a tensor 
product multiplicity is non- vanishing. The focus there is on gl(N) (and therefore also on su(N)) 
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and results are provided in the form of polyhedral combinatorial expressions. Our prescription 
described above is in general more explicit than results previously obtained. For lower rank 
su(N), though, the explicit inequalities may be obtained using various alternative approaches. 



4.1 A refinement 

Here we shall indicate how one may derive sets of inequalities determining when a tensor product 
multiplicity is greater than a given non-negative integer K 

?W > K . (33) 

The case K = was discussed above. 

Our approach is straightforward since the problem translates into studying when a convex 
polytope has a (discretised) volume bigger than K. To illustrate, let us consider su(3). In this 
case the volume is expressed as a single sum (p9|), so (|33|) is equivalent to requiring 



min{^i, A 2 , A 2 + - v 1 , A 1 + fi 1 - u 2 , -A 1 + A 2 + fi 1 - u l + v 2 , A 2 + // - fi 2 + u l - v 2 } 
- max{0, X 2 + fi 1 - fi 2 - u 1 , -A 1 + A 2 + /u 1 - v 2 } 

> K . (34) 



This leaves us the following 18 inequalities refining ( p0| ) 

K < Xi,fj,i,Ui , for i = 1,2 , 
max{A 1 -A 2 + ^ 1 -/i 2 ,-A 1 + / u 2 ,A 2 -// 1 } + i^ < v 1 < A 2 + fi 2 - K , 
max{-A 1 + A 2 -^ 1 +/i 2 ,A 1 -^ 2 ,-A 2 + ^ 1 } + K < v 2 < A 1 + /i 1 - K , 
max{-A 2 - /u 1 + /i 2 , -A 1 + A 2 - /i 2 } + K < v 1 - v 2 

< minjA 1 — [i 1 + fi 2 , 

-A 1 + A 2 + //} - K . (35) 

To the best of our knowledege, this is a new result. 

In the case su{4) the situation is already much more complicated. That is because the 
polytope is three-dimensional, and we cannot immediately use the triple-sum formula (|3l|). 
We recall that our simplifying method for obtaining an appropriate order of summation, may 
include redundant summations contributing zero to the final expression. We might therefore 
lose crucial information if we only considered the multiple sum formula. The remedy is to 
consider the original polytope, and require the defining faces to embrace a volume of at least 
the desired value. We would then be led to consider three-dimensional partitions of K+l, which 
is beyond the scope of the present work. For lower values of K the problem is straightforward, 
though. 



5 Conclusion 

By virtue of virtual BZ triangles we have obtained a polyhedral combinatorial expression for the 
su(N) tensor product multiplicities, different from the ones discussed in Refs. |], The main 
merit of our expression is that it admits a simple measurement of the convex polytope volume 
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in terms of a multiple sum formula. The latter is then a new and explicit way of expressing the 
tensor product multiplicities of su(N). 

As an application, one may derive explicit bounds on a triple of weights determining when 
the associated coupling to the singlet exists. To illustrate, the bounds were written for su(3) 
and sit(4). Also included was a brief discussion on how to generalise this to bounds describing 
> K, 

We believe that our multiple sum representation of the tensor product multiplicities provides 
a significant computational improvement over previous (combinatorial) results. In particular, 
it is expected to lead to considerable simplifications when implemented in computer programs. 

It is our hope that our results may find applications to the computation of fusion rules 
in conformal field theory with affine Lie group symmetry, the so-called WZW theories. Since 
tensor product multiplicities correspond to the infinite-level limit of fusion multiplicities, it is 
helpful to have simple descriptions of the former in order to understand the latter. 
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